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abstract The so-called STIT tessellations form the class of homogeneous (spa- 
tially stationary) tessellations of M. d which are stable under the nesting/iteration 
operation. In this paper, we establish the strong mixing property for these 
tessellations and give the optimal form of the rate of decay for the quantity 
\V(A n Y = 0, T h B n Y = 0) - n Y = 0)P(S n Y = 0)1 when A and B are 
two compact sets, h a vector of R d , Th the corresponding translation operator 
and Y a STIT Tessellation. 
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1 Introduction and notations 

Random tessellations, or mosaics, form an important class of objects of stochas- 
tic geometry. They have proven to be a useful tool for modeling geometrical 
structures appearing in biology, geology, medical sciences, and many others. 
Poisson hypcrplanc tessellation, Poisson Voronoi tessellation, and its dual tes- 
sellation Poisson Delaunay tessellation, are the most celebrated and tractable 
models investigated until now, and all are defined using a Poisson point process 
on an appropriate space. In the 1980's, Ambartzumian had the idea of an op- 
eration of iteration among the mosaics of M. d , namely the iteration, also called 
nesting operation. When regarded in R 2 , the nesting operation leaves stable the 
class of T-noded tessellations, while the superposition leaves stable the class of 
X-noded tessellations, such as the Poisson line tessellation. Nagel, Mecke and 
Weiss ([5], [6], [7]) have introduced the model of the STIT tessellat ion, following 
the idea of Cowan [3], which is a very natural object regarding this operation, 
because it possesses the feature of being STable under ITeration. This object 
can be used as a model for crack patterns, such as those seen on old pottery or 



*Univ. Lille 1, raphael.lachieze-rey@math.univ-lillel.fr 



1 



during earthquakes. 

The tessellation is homogeneous, i.e space stationary, and a proper choice of 
paramaters can make it isotropic, which yields a very interesting model. Many 
geometrical features of STIT tessellations have been investigated, including mo- 
ments of variables related to typical faces of the tessellation. Cowan, in [T] and 
[2], emphasized the importance of ergodic properties for random tessellations. 
In this paper we establish that all STIT tessellations possess the strong mixing 
property, which implies its ergodicity. Namely, if A and B are two Borel sets, 
and Y is the closed set of boundaries of the cells of a STIT tessclation, then 

P(Y n A = 0, Y n T h B = 0) - P(Y n A = 0)P(Y n B = 0) > 0, (1) 

| \h\ | — >oo 

where ThB is the set B translated according to the vector h. Moreover, when 
A and B are compact sets, we are able to give optimal rate of decay, of order 

PIT" 

After introducing the notations, we will give a brief decription of the construc- 
tion and the properties of STIT tessellations, and we will finally establish its 
mixing properties. The last section is devoted to the proofs. 

We introduce the following notations: 

• E = M. d is the d-dimcnsional cuclidcan space. 

• For any set f2 equipped with an unambiguous topology, B(Q) is the a- 
algcbra of Borel sets of SI. 

• T is the class of closed sets of E. 

• Q is the class of open sets of E. 

• K. is the class of compact sets of E. 

• H is the set of affine hypcrplanes of E. 

• l S (i ~ 1 is the unit sphere of E. 

• M. + is the set of all non-negative numbers. 
For A a subset of E: 

• int(^4) denotes the interior of A. 

• d(A) denotes the boundary of A. 

• conv(^4) denotes the convex hull of A. 

• span(^4) is the smallest linear subspace of E containing A. 
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For a scalar A > and a vector h £ K , the scaling and translation operators 
(resp. x — » A.t and T% : x — * x + h) are naturally lifted to J- , with the same 
notations. 

Since we are interested in hitting and missing probabilities here, we introduce 
the corresponding families. 

Definition 1. For A £ B(E), 

t a = {Cet I cnA^0}, 
F A = {Ce T I CnA = 0}. 

Definition 2. TTie toplology of closed convergence on T , or Fell topology, 
is the topology generated by the T K , K £ K. and the To, O E Q . 

For any set V C J 7 , we can define the Fell topology by restriction and 
the corresponding Borel cr-algebra B{V). We know that the c-algebra B(!F) is 
generated by the T K ,K £ K, (see [4]). Let us give definitions for tessellations 
and random tessellations. 

Definition 3. We define a tessellation of M. d as a countable set R of convex 
polytopes of M. d which satisfy 

(i) V C £ R, int{C) ± 0. 

(ii) U ceR c = m d . 

(iii) V C,C €R, C ^ C => int(C) f| int(C) = 0. 

(iv) V K &K., card(T K f| -R) < oo. 

Mc = {RE VQC) I R satisfies (i), (ii), (iii), (iv)}. 

Wo define the application y, from M r to J r , by 

cei? 

the closed set formed by the boundaries of all cells of the tessellation. 

Wc can define the reciprocal application: For Y a closed set, we define TZ(Y) 
as the set of the topological closures of all connected components of Y's com- 
plementary. 

Definition 4. Let M. be the class of tesselations as closed sets ofR d : 
M = {Y E T | K(Y) E M c }- 
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Let B(A4 C ) be the cr-algebra on A4 C that makes 

y : (M,B(M)) -f (M C ,B(M C )) 

measurable (where M. is endowed with the restriction of the Fell topology). 
Since the applications y and 1Z are measurable reciprocal bijections, we can 
define a random tessellation as a random variable taking values in either Aic or 
M, and henceforth we will often jump from one description to another, without 
bothering with the measurability questions. 

According to Choquet's theorem (see [3], [5]), the law of a random closed set is 
characterized by its capacity functional, which is defined by: 

Definition 5. Let Y be a random closed set of law P on B(J-). We define 



V K e /C, T Y (K) = P(F K ). 



If one is given two random closed sets Y and Z, one will be able to say that 
they have the same law if there is equality of their capacity functionals. 

y = Z«VA'eK, Ty(A') = T Z (K). 

2 STIT Tessellations 

We will first give the decription of the local construction of a STIT tessellation, 
and then use spatial consistency to extend the definition to the whole space. 
We will then define the nesting iteration operation and give its connections with 
STIT tessellations, and finally give the main properties of such tessellations. 
A STIT tessellation, like the Poisson hyperplane tessellation, is defined from 
an hyperplane Poisson process. All details oncerning the construction of STIT 
tessellations can be found in [5], [5] and [7]- 
Let us begin with some definitions. 

Definition 6. For A 6 B, we set 

[A] = nn t a 

the hyperplanes hitting A. 

For B another borelian set, we also define [A\B] as the set of all hyperplanes 7 
that separate strictly A and B, i.e such that A and B are contained in different 
open half- spaces of E\j. 
We write in short, for x, y € M. d , 

[x] := [{x}}, 
and [x\y] := [{x}\{y}]. 
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We equip Ti. with the Borel a-algebra induced by Fell topology and let A be a 
translation invariant locally finite measure on Ti. 

Given 7 £ Ti \ [0], we denote j + the closed half-space delimited by 7 not con- 
taining 0, and 7~ the other one. Note that A([0]) = since A is stationary and 
locally finite. 



Let a be a positive number. 
We identify Ti with R + x S 1 where an element (r, u) of M + x stands for 
the hypcrplanc 7 at distance r of the origin and normal exterior vector u (i.c 
normal vector directed towards 7 + ). 

Since A is assumed to be stationary, we can write, along with this identification 

A = A + ® v, 

where A + is the Lebesgue measure restricted to R + and v is a finite measure on 
We assume that v has the following property: 



span(supp(z^)) = E. 



(2) 



The latter property ensures that tessellations generated by the hyperplane pro- 
cesses with distribution A will a.s have compact cells. 

On a compact window W £ K, with non-empty interior, we are going to define 
the tessellation with parameters a, A as a random time process taking values in 
M, with a the time parameter. Let B be Young's binary infinite tree, where 
the leaves are labeled so that the top leaf is attributed to number 1, and the 
fc-labeled leaf has the leaves 2k and 2k + 1 as daughter leaves. We now attach 
to each leave a pair (efe, 7&), where (efe)fe>i is a family of iid random exponential 
variables with parameter A([W / ]), and (7fe)fe>i are iid random hypcrplanes with 

laW A[W] ■ 

The tessellation is defined as a process of cell division, where each cell is iden- 
tified with a leaf of the tree, and the first cell is the window C\ = W itself. 
We describe it in term of a birth and death process, where a cell dies when it 
is divided in two daughter cells. The death times (dk), birth times (bk), and 
daughter-cells of (Ck)>x are defined recursively in the following way ([s\ is the 
integer part of s): 

61 =0 
bk = dik/2\ 
dk = b k + £fc 

c 2fe = c k n 7- 
c 2 fe+i = c k n 7^ 



We then define 



R a ,w 

y a ,w 



{C k I b k <a,C k 7^0}, 
|J dC k . 



b k <a 
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We call in the sequel P a ,w the law of this tessellation. Note that A is an implicit 
parameter of the model. 
We also write in short 

^V»,w := T a ,Wi U a ,w = 1 — T a< w 

Remarks 

1. The only non-intuitive feature of this construction is that 7*. might not 
be in [Cfc], which would make one daughter cell be and the other be 
Ck itself. There is no theoretical objection, but one can remedy this by, 
instead of choosing an iid family (efe), attach independantly to each cell 
Ck the death rate AQCfc]) and an hyperplane drawn from ^rgj-fp , so that 
the law is not modified and each hyperplane indeed hits the cell to which 
it is attached. 

2. When A is isotropic, the death rate of a cell C, A([C]), is proportionnal 
to the perimeter of C. 

3. We mentioned in the introduction a construction from a Poisson hyper- 
plane process. Let 6 be a Poisson process on R + x TC with intensity 
H[W])X+ ® A( |$p . Since a.s, for all t > 0, card(({t} x H) n 6) < 1, 
we can define a random sequence (Tk,jk)k>i such that, almost surely, for 
all k > l,T k < Tfc+i, 9 = {(r fe ,7 fe ) | k > 1} and (r k+1 -Tk,Jk) satisfy the 
hpothescs of the previous construction. 

Consistency 

Nagel and Weiss [7j have established the consistency property: 
Theorem 1. If W C W are two compact sets with non-empty interiors, then 

Y a . w = {Ya..w n^)U dW. 

In particular, for any compact set K C int(M / ), 

T a , w {K)=T a , w ,[K), 

and we can define on K, 

T a {K): = T a . w (K) 

U a {K) : = U a w{K) for any W containing K in its interior. 

Due to Choquet's Theorem, this allows us to define a tessellation Y a on the 
whole plane, with capacity functional T ai such that 

WW &IC, (Y a n W) U dW C = ] Y a , w , 

which is called the STIT tessellation with parameters a and A. We call P a its 
law on M. or, equivalently, on M c - Mecke, Nagel and Weiss [6] also provided 
an explicit global construction on E. 
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Iteration 

Let us now define briefly the nesting operation. A more elegant definition, 
with the help of marked hypcrplane Poisson processes, is given in [5]. Given two 
random tessellations Y and Y' having laws Q and Q 1 , we define the law Q EE Q' 
as that of the tessellation Z obtained in the following way: 
Let us take a copy of Y and an independant iid family {Y q ) q >\ of copies of 
Y'. Let us put on the cells of Y a measurable labelling independant of the Y£, 
Y = Uq>i 9C q , and let us define Z as the random closed set: 



so that Z is obtained by nesting an independant copy of Y 1 in each cell of Y. 
WenotcZ = y®(y 9 ') 9 >i. 

A random tessellation is said to be stable under iteration iff its law P satisfies 



It has been shown in [7J that STIT tessellations are the only random tessellations 
stable under iteration. Moreover, the class of STIT tessellations is attractive for 
the iteration in the following sense. Let us take a random stationary tessellation 
Z with law Q, such that a is the meanlength of Z and A is its directional 
distribution. We define the sequence of laws of tessellations (Q n )n>i by: 



Hence, the law Q n is that of Z iterated n times with itself. It is shown in 
[7J that the normalized probability Q n {^-) converges weakly to the law of the 
STIT tessellation with parameters a and A. So, the STIT tessellation is really 
the central object regarding this nesting operation. In a way it is not surprising 
since it is itself constructed in each compact window as some kind of iteration 
of the tessellation generated by a single hyperplane. 



3 Strong mixing property 

In this part we will prove relation |T]). and give a rate of decay close to the 
optimal for compact sets. 

Here A = A + ® v is a measure on TL satysfyings assumption and a is a 
parameter strictly positive. 

In this section the probability space is (A4, B(A4), P a ), and Y a is a random 
variable with law P a . 




p 





Q, 

Qn EB Qi- 
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Consider the set T = {T/j | h G M. d } of all translations, seen as operators on 
T . The stationarity of a random closed set X with law P ensures that for all 
T G T, P is invariant under T. Now we call T-invariant set of B(J-) every Borel 
set C of the cr-algcbra B(J-) such that 

V T G T, TC = C. 

For instance /C, 7i, or "the class of all tessellations having a cube as one of 
their cells" are invariant sets. Given a stationary law P on T , the dynamical sys- 
tem (T ', £>(T), P, T) is said to be ergodic if every T-invariant set has probability 
or 1. and strongly mixing if for all C, C in J-, 

P(C,T h C) heE ) P(C)P{C). (3) 

| \h\ | — >oo 

When X is a random closed set with law P, we simply speak of the stationarity 
the ergodicity or the mixing property for X. The mixing property ensures the 
asymptotic independance of relatively far away sets. Since any T-invariant set 
C satisfies P(C,ThC) = P(C), the strong mixing property implies ergodicity. 
According to Lemma 9.3.1 in [8], it suffices to show ([3]) for sets C,C drawn 
from a semi- algebra A generating C. That is why, in order to obtain the strong 
mixing property for Y, it is sufficient to show the following theorem: 

Theorem 2. For all A, B e B, 

P a (T A nT T - B ) P a (T A ) P a (T B ) . 

1 1 h 1 1 — >■ oo 

Note that T h T B = T T ~ %B ' , but we state the theorem with h instead of — h 
for more simplicity since the role of h is symmetric. Since the {J- k )kek a ls° 
generate the Borel cr-algebra of the Fell topology, still using Lemma 9.3.1 in [8], 
it suffices to show it for A,BelC, and the proof is given in section [H We are 
able to indicate the optimal rate of decay for that we first need to introduce a 
function. 

Proposition 1. Given u G S^ 1 , the measure of all hyperplanes separating 
and u has the following representation: 

C(u):=A([0|u])=A([0,u]) = i / v(dv)\{u,v)\, 

£ Js d - 1 

and so 

Wu,veS d -\ \({u)-((v)\<^\\u-v\\. 
Furthermore, there exists k > such that 

VuG 5 d_1 ,C(«) > K - 
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This result, shown in section 0] along with others, is meaningful regarding 
the asymptotic behaviour of expression {T]), given in the following theorem: 

Theorem 3. For all A, B £ /C, There exists a constant x(A, B, a) > 0, invari- 
ant to rigid motions of A and B, such that, 

|P.^n^)-P.(^P.(^l < (4) 

Furthermore, when A and B each have a single connected component, 
P a {T A ,T T - B ) 



Pa(T A )Pa{T B ) 



1 



IWIC(VIHI) 



+ (5) 



The asymptotic dependence of T A and T ThB is closely connected to the 
probability that the cell W = conv(AUThB) is hit before time a by a hypcrplane 
that separates A and T^B, which dependency is asymptotically in | \h\ \((h/\\h\\), 
using the stationarity of A (see lemma [3]). 

According to ([5]), ^ is optimal in the sense that the inequality is an equality 
when A and B are connected sets. When they are not, we think that we are close 
to optimality, but the estimation is more difficult because the capacity functional 
has a much more complicated expression, even with only two components (see 
[7]). To deal with it, we have to estimate the lipschitziannity of the capacity 
functional as a temporal function of a. We obtain the following expression, 
worth to be noted: 

Let t > be a real positive number, and K a compact set. Then 

There exists A A > > 0, < T a+t (K) - T a {K) < t\ K . a . (6) 

Furthermore, A A ,a is invariant to rigid motions of K . 
If one had a non-zero linear lower bound in ([S]), one could probably, by a slight 
modification of the proof of theorem [3J obtain a lower bound of the left-hand 
term in ((4]) by one of the form of the right-hand term. 

4 Proofs 

Proof of © 

Lemma 1. Let < t, a be real positive numbers, K a compact set, and 
X K , a = A([conv(K)])(l + aA([conv(K)]))(l - T a(K) ). Then 

0<T a+t (K)-T a {K)<t\ K ,a- 



We use the notations of section ® with W = conv(K). Let us call N a the 
number of times conv(if ) has been hit by a hyperplane up to time a: 

A^ a = card{ k | 7fc G [C k ]}. 
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For n G N we set 7r„(o) = P(A a = n) > 0. Almost surely, conditionaly to 
(N a = n), R a _w contains n + 1 cells Ci, C„+i with non-empty interior. We 
call A^ := K n Q. We have 

T a+t (K) - T a (K) = U a {K) - U a+t (K) = U a (K)F(K n Y a+t ± 0|A R F a = 0) 
= £/ a (A) ^ 7r n (a)P(3i, 1 < i < n + 1, tf< n F a+t ^ 0|A a = n, if n Y a = 0) 

n>0 

n+1 



< u a {K) 7r n (o) ^ W n r a+t ^ 0|A a - n, if n r a = 0) 

n>0 i=l 
n+1 

C/ a (A) ^ 7T„(a) ^ P(conv(A l ) n Y a+t ^ 0|A a = n, conv(A,) n Ya 



n>0 i=l 

At this point, for conv(Ai), knowing that Ki is in the interior of a cell at time 
a amounts to know that conv(Ai) has not been hit up to time a, and so, due to 
the lack of memory and the consistency property, 

P(conv(A 4 ) n Y a+t ± | N a = n, conv(A 4 ) n Y a = 0) = 1 - [/ t (conv(A l )) 

< 1 - £/ t (conv(A)) 

= e -tA([conv(KT)])_ 

This yields 

n+1 



U a (K)-U a+t (K) < U a {K)Y,Kn{a)Y,[ l 



e -tA([conv(if)]) , 



ti>0 i=l 
n+1 



n>0 i=l 

< U a (K) »r„(o)(n + l)iA([conv(A)]) 

Tl>0 

< [/ a (A)iA([conv(A)])(l + E(JV )). 

The number of hypcrplancs involved in the cell division process with initial cell 
conv(A) up to time a is less or equal to a Poisson variable with parameter 
aA([conv(A)]) (remember that some hypcrplancs fall outside the cell they were 
supposed to divide). In consequence we have 

E(N a ) < aA([conv(A)]), 

which concludes the proof. 

□ 
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Proof of Proposition Q] 

For u £ iS d_1 we have by stationarity 

((u) = A(]0,u[) = A([0,u]) = -A([-u,u\) by stationarity. 

Now, given a hypcrplane 7 at distance r from the origin with exterior normal 
vector v, 

7 6 [— u, u] 4=> < r, 

and so 

^ JS d - 1 

Yovu^u' ^ S d -\\C{u) - au')\ < If I \(u,v)\-\(u',v}\ \u(dv) 

< - \(u,v) - (u',v)\v(dv) 

\\u-u'\\ f 



d - 1 



\u-u\\ d _ u 



Let u £ Thanks to assumption ^ is not concentrated on u^, hence 

J S d-i \ ( u , v)\v(du) > 0. Since is compact, we have the existence of k > 

such that C(u) > k. 

□ 



Proof of Theorem [3] 



We are first going to establish some inequalities for A, B any compact sets, 
and then add a drift h to give an upper bound when expressions become too 
complicated. 

Let W = conv(A U B). 
The key is to consider the tessellation valued time process defined in section [21 
(Xt,w)t>o- We keep the notations e\ = inf{i | Y 0t w 7^ and 71 is the first 

hyperplane dividing C\ = W. 
We introduce the event 



Fa.b = (71 6 [A\B],E! <a). 

When A and B are far away from each other, Ta,b is very likely to happen, and 
after Ya,b occurs the behaviours of A and B are independants because they are 
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in disjoint cells. 

We first have to estimate 



\p a (F«,F*,r AtB )-p a (F«,f»)\ < p a {v% B ) 

< e -ah([W])_ 

In consequence, we can show formulae (gj| and © with P a {T A , T ThB , T A . ThB ) 
instead of P a (T A , J rThB ), because their difference is negligible with respect to 
the expected decay rate. Rigorously this gives: 

P a (T A ,T B SA. B ) = / P a (ei G di, 7l G [A\B],T A ,T B ) 
Jo 

= I Pai^i G di,7i G [A\B])P a (J 7A , T B \e\ =t,ji G 
Jo 

When e\—t and 71 G [j4|P], it means that A and B have not been hit up to time 
t and are both still contained in a cell (but the cell encapsulating A is different 
from that of B). Thanks to the consistency property, and the independent 
behaviour of distinct cells after their birth, we have 

P a {T A ,T B \d = t, 7 i G [A\B]) = U a -t(A)U a - t (B). 
Since the sequence (7^) is independent of (e^), we have 

P a (ei Gdi.7! G [A\B]) = Paid G di)P a ( 7l G [A\B]) 

U Jj A([W]) ' 

and one finally has 

P a (^ A ,^ s ,r AB ) = A([A|B]) r e- tA U w »U a - t (A)U a . t (B)dt, (7) 

Jo 

so that 

I P a [T A , .P s , T A , B ) - P a (T A )P a [T B ) I 

= A([A|P]) /V' A (^% Q _ t (,4)[/ a _ t (P)di- 

/>oc 

U a (A)U a (B) / e-* A ([ w ])A([^])dt 



A([A|P]) / e -*A(W> |C/ a _ t (A)[/ Q _ t (P) - [/ a (A)t/ a (P)| dt+ 
e - iA a w ])c/ a (^)C/ Q (B)A([W])dt + 

/"OG 

+ A([A|P]) / e-^^Mi. (8) 

</ a 



A([A|B]) x 



A([W]) 
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Now it is time to consider the drift of B, in order to make a proper upper bound 
estimate. We will use T^B instead of B in the previous formulae and use the 
following facts. 

Lemma 2. Let A,B £ K, and h be an element o/R d . We set 
W h = com(A,T h B). Then 

IWI i 

e- u du = o(_), (9) 

A([^]) = A([A|T,S])(1 + (1)), (10) 

A([w„]) = v (VIHI)IHI(i + o(i)). (ii) 

Proof. Let us prove the two latter formulas: Let a € A, /3 e B. If 7 is a 
hyperplane separating a and (3 — h and hitting neither conv(^4) nor conv(T^i?), 
it also separates A and T^i?. Hence 

\A([A\T h B}) - A([a,(3~h})\ < A([conv(A)]) + A([conv(T, l B)]) 

< A([conv(A)]) + A([conv(B)] = o(\\h\\). 

For similar reasons, 

< A(W h ) - A([A|T h B) < A([conv(A)]) + A([conv(B)]) = o(||ft||), 
and we only need to prove 

A([a, /?-/*]) = ^(VINDINKi + oa)), 

which is given by the following result. 

Lemma 3. If fi — a = — a\ \u, with u £ , 

A(M/3]) = - aHC(«)- 
Proof. It suffices to show that for all e > 0, u € S x , n G N, 



A([0|(ra + l)eu]) = A([0|neu]) + A([0|eu]), 

because, since A is stationnary and locally hnite, V x £ R d , A([x]) = and then 

we will be able to obtain the result by induction and approximation. 

Since [0|(n + l)eu] is the disjoint union of [0|neu], [neu\(n + l)eu] and [new], we 

have 

A([0|(n + l)eu]) = A([0|neu]) + A([neu\(n + l)eu\) 
and the result follows by stationnarity of A. 

□ 
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Finally we can end the proof of the lemma: 

A([W h }) = \\a-P + h||C( ,, a ~flta,| )(l + 0(1)) = ||ft||C(VII^II)(l + 

\\0i-p + h\\ 

since £ is continuous. 

□ 

Let us define 

5(A,B,a) = \A,a + ^B,a- 

> >>A,a-t + >^B,a~t for all t > 0. 

We set 

x(A,B,a) = 8{A,B,a) + A([conv(A)]) + A([conv(P)]). 
Using ((6|) and |(5J), and the previous lemma, we get 

\p a (F A ,r ThB ,r A , ThB ) - p a (T A )p a (T B )\ 

< A.([A\T h B])S(A, B, a) [ e~^m + (1 - ^j^) + o(^) 

g(A,B, a )A([A|r h B]) A([conv(A)])+A([conv(T^)]) 1 
A([^])2 + A([W fc ]) Vir 

x(A,P,a) 

^(VINI)INI (1 + 0(1)) ' 

and formula ^ is proved. 

When ^4 and B are connected, thanks to the simple expression (see [7J) 

U u (A)=e- uA ^\ 
using ([7]) we have the exact expression 

\p a (T A , t t - b , r A , ThS ) - p a (^ A )p a I 

= P a (T A )P a (T B ) p e -KMW h )-A(A)-A(B)) dt 

oo 







Hence 

p a (^,^ B ,r^ ThB ) 



- i 



\hMh/\\h\\) 



(l + o(l)), (12) 



P a (F A )P a (F B ) 
and (|5|) is proved. 

□ 
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